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Abstract 

We consider a simplified physics of the could interface where conden¬ 
sation, evaporation and radiation are neglected and momentum, thermal 
energy and water vapor transport is represented in terms of the Boussi- 
nesq model coupled to a passive scalar transport equation for the vapor. 
The interface is modeled as a layer separating two isotropic turbulent re¬ 
gions with different kinetic energy and vapor concentration. In particular, 
we focus on the small scale part of the inertial range as well as on the dis¬ 
sipative range of scales which are important to the micro-physics of warm 
clouds. We have numerically investigated stably stratified interfaces by 
locally perturbing at an initial instant the standard temperature lapse 
rate at the cloud interface and then observing the temporal evolution of 
the system. When the buoyancy term becomes of the same order of the 
inertial one, we observe a spatial redistribution of the kinetic energy which 
produce a concomitant pit of kinetic energy within the mixing layer. In 
this situation, the mixing layer contains two interfacial regions with op¬ 
posite kinetic energy gradient, which in turn produces two intermittent 
sublayers in the velocity fluctuations field. This changes the structure of 
the field with respect to the corresponding non-stratified shearless mix¬ 
ing: the communication between the two turbulent region is weak, and 
the growth of the mixing layer stops. These results are discussed with 
respect to experimental results with and without stratification. 


1 Introduction 

Warm clouds as stratocumuli swathe a significant part of earth’s surface and 
play a major role in the global dynamics of atmosphere by strongly reflecting 
incoming solar radiation - thus contributing to the Earth’s albedo - so that an 
accurate representation of their dynamics is important in large-scale analyses 
of atmoshperic flows [T] . They are controlled by the tight interplay between ra¬ 
diative driving, turbulence, surface fluxes, latent heat release, and entrainment. 
Among them, the mixing and entrainment processes at the cloud top have been 
identified as fundamental to determine the internal structure of warm clouds, 
so that a clear and complete understanding of their physics is required [5]. As 
pointed out by Malinowsky et al. [3] , data from most field campaigns and large- 
eddy simulations are too poorly resolved to allow to infer the details of the 
interfacial layer, even if they indicate that, in order to allow for entrainment, a 
high level of turbulence must be present. For this reason, in this work we study 
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Figure 1: Kinetic energy spectra. 
Contextualization of present study 
(black spectra, small inertial and 
dissipative range) respect to spectra 
from in-situ atmospheric measure¬ 
ments [311 EH] (colored spec¬ 
tra, energy injection and low wave- 
number inertial scales). 


the local transport through a clear air/cloud interface through DNS (Direct 
Numerical Simulation). As our focus is on the dynamics of the smallest scales 
of the flow which influence the microphysics of warm clouds, we have simulated 
an idealized configuration to understand, under controlled conditions, some of 
the basic phenomena which occur at the cloud interface over length scales of 
the order of few meters. In these conditions, we solve scales from few meters 
down to few millimeters, that is, we resolve only the small scale part of the in¬ 
ertial range and the dissipative range of the power spectrum in a small portion 
(6 TO X 6 TO X 12 to) of the atmosphere across the clou - clear air interface. This 
allows us to investigate the dynamics of entrainment which occurs in a thin 
layer at the cloud top, which a smaller scale with respect ot the scale explicitly 
resolved in large-eddy simulations of clouds [4]. In this preliminary work, we 
focus on two concomitant aspects of the cloud top mixing layer: the effect of 
the presence of a stratification and of a turbulent kinetic energy gradient. We 
do not consider the wind shear neither the phenomena linked to the processes of 
evaporation and condensation and radiative cooling which are important in con¬ 
ditions of buoyancy reversal [3E]. Therefore, our simulations were performed 
by applying the Boussinesq approximation to the Navier-Stokes momentum and 
energy equations together with a passive scalar transport equation which models 
the water vapor transport. 

The details of the physical problem we have considered and of the governing 
equations are given in section 2. Section 3 contains a selection of our main 
results about intermittency, energy redistribution and entrainment. Conclusion 
remarks are in Section 4. 


2 The physical problem 

We consider the interface between clear air and cloud ina6TOx6TOXI2TO 
parallelepipedic domain. Compared to in situ measurement of the atmospheric 
energy spectra, as shown in figure [1] we are able to simulate the lowest part 
of the inertial range and the dissipative one. As shown in figure the system 
is composed by two homogeneous and isotropic turbulent regions that interact 
through a mixing layer, whose initial thickness has been set of the same order of 
the integral scale of the turbulence background here assumed equal to 3 • 10“^ 
m. The two isotropic regions have a different kinetic energy and we assumed 
that the kinetic energy is higher in the cloud than in the external region. The 
root mean square of the velocity inside the cloud is Urms = 0.2 m/s, and the 
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Conipntatioiial domain: ~ 12 m 

Initial cloud/clean air interface: Ag ~ 0.3 m 
Computational grid: 1024^ x 2048 


standard temperature lapse rate 




Top interface 


E 2 - Lower kinetic energy 


unstable 


El - Higher kinetic energy 


stable 


CUMULUS 2 r; 1000 m 
water vapor concentration represented 
as a transported passive scalar 


Figure 2: Scheme of the initial conditions. Ei is the mean initial turbulent 
kinetic energy of the bottom region (inside the cloud), E 2 of the upper re¬ 
gion (outside the cloud). For the top cloud mixing here presented we consider 
E 1 /E 2 = 6.7. The water vapor is initially present only inside the cloud (bot¬ 
tom region). The zoom in the yellow circle is an example of initial temperature 
perturbation 6' of the standard boundary layer lapse rate. 

energy ratio between the cloud energy Ei and the external region energy E 2 
is equal to 6.7. This energy ratio is of the same order of the ones measured 
in warm clouds (see, e.g., [3]) and, furthermore, it allows us to compare our 
results with experiments on shearless mixing (see m 112] ') in absence of any 
stratification. In our simulations the Taylor microscale Reynolds number Re\ 
of the higher energy region is equal to 250. Buoyancy is taken into account 
through a local perturbation 9' in the profile of temperature distribution 9 
inside the troposphere, and the water vapor concentration y is considered as a 
passive scalar. The Prandtl (Pr= 0.74) and the Schmidt number {Sc = 0.61) 
considered refers to an altitude of 1000 m s.l. 

The initial conditions for the temperature perturbation is described in [5] and 
in table [TJ The ratio between inertial and buoyancy forces is expressed by the 
Froude number Fr, based on the maximum gradient within the initial interface, 
which ranges from 31.2 (negligible stratification) to 0.62 (strong stratiheation). 

We use the continuity, momentum and energy balance equations within the 
Boussinesq approximation, which holds for small temperature variations |13| : 


V • u' = 0 

Ti 

— -1- (u' • V) u' = —V- -1- -I- ag0' 

at p 


( 1 ) 


( 2 ) 
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Table 1: Initial stratification level parameters. G is the maximum gradient of 0, 
expressed in terms of the standard troposphere lapse rate Go = 0.0065 Km~^\ 
is the characteristic Brunt-Vaisala frequency of initial con- 
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dition. The Froude number Fr= and the Reynolds Buoyancy number 

Refe = give a measure of the order of magnitude of the buoyancy forces 
compared with the inertial terms (e is the initial energy dissipation rate, v the 
kinematic viscosity). 
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where the temperature 9 is composed as the sum of a fluctuation d'(x, t), 
a static component 6 {x^) = Gxs, and a reference constant temperature 9q, 
p = p + agx 3 (do + GX 3 / 2 ) is the total hydrodynamic pressure (where p is the 
fluid-dynamic pressure, a the thermal expansion coefficient, g the gravity accel¬ 
eration), u' is the velocity fluctuation and x is the vapor concentration of the air 
- water vapor mixture present in the cloud, here considered as a passive scalar. 
The constant k and are respectively the thermal and water vapor diffusivity. 

The simulations are performed using our home produced computational code 
that implements a pseudospectral Fourier-Galerkin spatial discretization and an 
explicit low storage fourth order Runge-Kutta time integration scheme. Eval¬ 
uation of non-linear (advective) terms is performed through the 3/2 de-aliased 
method M- The initial conditions for the velocity field are obtained by a linear 
matching of two different isotropic homogeneous turbulent fields (that are ran¬ 
domly generated, respecting physical conditions imposing spectra, solenoidality, 
integral scale and kinetic energy) [15]. The grid has 1024 x 1024 x 2048 points, 
and allows to capture all the turbulent scales from the greatest (integral scale i) 
to the smaller (Kolmogorov scale rj). The computational code uses a distributed 
memory paradigm through the MPI libraries: the simulation were performed 
at the TGCC Curie supercomputer within the PRACE project n° RA07732011 
for a total of 3 million cpu-hours. 


3 Results 

According to the ratio between buoyancy force and kinematic forces (that are 
advection and diffusion), the evolution of the system can be split in two main 
stages. As long as the ratio remains small, there are no significant differences 
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Figure 3: Spatial distribution along the 
vertical coordinate of the kinetic en¬ 
ergy (a), temperature variance (b) and 
water vapor concentration variance (c) 
at t/r Ki Q for different levels of strati¬ 
fication. 


with respect to a non-stratified case. On the contrary, as the stratification per¬ 
turbation level become higher, buoyancy effects are no more negligible: differ¬ 
ences are present from both a quantitative and qualitative point of view. These 
considerations can be observed through the statistical analysis of the simulated 
fields. The statistics are computed by averaging the variables in the planes nor¬ 
mal to the mixing direction (with a sample of 1024 x 1024 data-points), focusing 
on the variation along the vertical (non-homogeneous) direction. The effects of 
the different stratification levels are clearly visible on the second order moment 
of velocity, temperature, and vapor as shown in figure |3] (a - b - c). When the 
stratification level is mild {Fr > 4) there are no relevant differences with re¬ 
spect the neutrally buoyant flow, while significant differences appears for intense 
stratihcation {Fr < 2). In particular, in correspondence of the local tempera¬ 
ture perturbation, the formation of a layer with a pit of kinetic energy can be 
observed. The presence of such a layer deeply changes the physics of the sys¬ 
tem, because in this situation two interfaces are produced. The first interface, 
(which is present also in the absence of stratification), now separates the high 
turbulent energy region from the pit, while the second one (not present without 
stratification) separates the low turbulent energy region from the center of the 
mixing layer. Therefore, a strong stratification induces a physical separation 
between the two external regions, greatly decreasing the interaction between 
them. Both interfaces present an intermittent behavior, as shown in figure]^ (a 
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Figure 4: Spatial distribution in the vertical direction of the skewness and 
the kurtosis of vertical component of the velocity (a - b), temperature (c - d) and 
water vapour concentration (e - f) at t/r = 6 for different levels of stratification. 
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- b) by skewness and kurtosis distribution (respectively third and fourth mo¬ 
ments normalized with the local variance). In fact, two peaks of skewness and 
kurtosis can be observed in the highly stratified case: one is placed inside the 
cloud, and the other is placed close to the position of the intermittency peak in 
case of absence of stratification, see data at Fr = 0.6 in figure]^ (a). Observing 
the magnitude of the kurtosis maximum in figure |^(b), it can be noted that the 
peak inside the cloud reaches values as high as 4, that are about the 10% larger 
than when the stratification is milder. 

The behavior of the temperature and the water vapor fields appear qualita¬ 
tively analogue as far as the shape of the statistical distributions are concerned, 
but are quantitatively different. See, for instance, the temperature and the 
vapor concentration peaks in the skewness and kurtosis distribution, shown in 
figure]^ (c - d - e - f). 

It can be also observed that higher levels of stratification produce a relevant 
reduction of intermittency in the flow, with a drop of about 70% in the peaks 
of skewness and kurtosis. The interaction between the two regions aside the 
interface is greatly reduced, so the fluctuation at the sides of the mixing layer 
are damped, preventing the formation of the intermittent layer typical of the 
passive scalar transport [MHz]- This strong reduction in skewness and kurtosis 
is coupled with a slight increase in the higher order moments of vertical velocity 

- for sufficiently strong stratification - which is in fair agreement with the trend 
observed in [H p. 158]. Moreover, observing both figures and it is clear 
that the thickness of the mixing layer is reduced in case of intense stratification 
(see next section for more details). 

3.1 The onset of a kinetic energy pit 

As shown in figure (a), in case of high stratification level, in the center of the 
domain - where the initial temperature step is placed - the onset of a layer with 
a kinetic energy lower than both the external regions can be seen. This layer can 
be considered as a pit of kinetic energy. Varying the stratification intensity, the 
genesis and the evolution of such pit can be measured by considering the ratio C, 
between the variance of vertical velocity in the low energy region u^mle (mean 
value) and in the center of the pit u^min (where the variance reaches its mini¬ 
mum); the temporal evolution of C is shown in figure [^(b). Qualitatively sim¬ 
ilar results have been observed in the large-eddy simulations of stratocumulos- 
topped planetary boundary layer carried out by several physics of the atmo¬ 
sphere research group, as described in [B fig.7(c) and 8(a) at pag. II]: in par¬ 
ticular, in the case of sufficiently strong stratification, the trend of our vertical 
velocity variance, in figure [51]a) , is analogue to those observed in the LES car¬ 
ried out by the NCAR group (Deardoff TKE model [50]) and the WVU (ARAP 
TKE model [H]). In this simulations they consider the planetary boundary 
layer with Rex ~ 5500 and Fr « 0.4. 

A visualization of such phenomenology is represented in figure where the 
vertical velocity fluctuations in a vertical slice of the domain are represented 
using an elevation plot (where such elevation is proportional to the square of 
rta). In presence of a mild stratification, Fr = 12, even after 8 time scale, there 
is a smooth mixing layer between the high (left) and the low (rigth) energy 
regions. The differences in the case of strong stratification, with Fr = 0.62, 
after 4 time-scales (b) and 8 time scales (c), are clearly visible: a separation 
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(a) 


(b) 


Figure 5: (a) Temporal evolution of the spatial distribution along the vertical 

direction 0:3 of the turbulent kinetic energy when Fr=1.8. A pit of energy 
appears after about 6 eddy turnover times in correspondence of the original 
interface, (b) Temporal evolution of the ratio between the mean vertical 
velocity variance in the lower energy (clear air side) region uz,mle and the 
minimum value of the vertical velocity u^^min- When this ratio departs from 
1 , a pit of kinetic velocity appears as shown in (a). 




(c) 


Figure 6 : Visualization of the vertical 
component of the velocity in a vertical 
plane. The elevation is proportional to 
the square of vertical velocity fluctua¬ 
tions and the colors respect the velocity 
directions (blue for downward , red for 
upward) in various case: Fr = 12 af¬ 
ter 8 eddy turnover times (a) and for 
Fr=0.62 after 4 and 8 eddy turnover 
times (panels b, c respectively). In the 
last two panels it can be observed the 
formation of the pit of kinetic energy. 











(a) (b) 

Figure 7: Thickness of the mixing layer of velocity (a) and temperature (b). For 
both fields, the mixing layer thickness is the same as defined in m- 


layer is present in the center of the domain, that becomes even more evident as 
the time pass by. 

As said, the presence of the pit generates a physical separation between 
the two external regions, by damping the turbulent mixing, and thus reducing 
the exchange of information. As a consequence, there is a saturation of the 
thickening of the mixing layer A^; such interruption of the growth is represented 
in figure 0 (a) 

Looking to the temperature mixing layer thickness Ag, shown in figure]^ (b), 
it can be seen that, for strong stratification, the thickening stops approximately 
after the same amount of times scale required by A^;. In that case, contrary 
to what observed for the kinetic energy, the thickening does not stop suddenly, 
but rather with a transient that lasts a couple of time scales. 

3.2 Entrainment 

The entrainment of clear air inside the cloud is an important aspect of the top 
cloud interface as it concurs in the evaporation/condensation of droplets inside 
a cloud [1]. 

In any plane parallel to the interface, in absence of a mean velocity, only 
downward velocity fluctuations can transport clear air into the cloud. Their 
presence can be represented by a marker function ij) that is equal to I where 
Us is negative, and 0 otherwise. Its average in each horizontal plane, shown in 
figure 1^ (a), shows a small deviation from the mean value of 0.5 which would be 
observed in an homogeneous and isotropic flow. This implies that upward and 
downward fluctuations are almost equally distributed; the maximum departure 
from an homogeneous distribution is about 4%, and the spatial distribution of 
Ip looks like the one which has been observed in the third order moment of the 
velocity, see figure (a). 

The entrainment of clear air is responsible of the growth of the cloud. In 
fact, the velocity = dxs^i/dt of the cloud top interface {{xs^i) is the mean 
vertical position of the cloud top, defined as the position where the mean va¬ 
por concentration is 25%) has often been used as a parameter to measure the 
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(a) (b) 



(c) 


Figure 8: Evolution of the entrain¬ 
ment across the top cloud. Fraction ^ 
of downwards velocity (a) and vertical 
variation of the mean flux of clear air 
into the cloud (c) when with Fr= 1.8 
Evolution of entrainment velocity 
normalized with the high kinetic energy 
El root mean square (b), see hgure 


entrainment rate 011]. In figure]^ (b), the temporal decay of for different 
levels of stratification is represented. In presence of weak stratification, that is 
Fr larger than 4, its value gradually decreases with an almost exponential law, 
due to the decay of the kinetic energy. On the contrary, when the stratification 
is stronger, the damping of We is much faster, and the entrainment vanishes af¬ 
ter few times scale, when the presence of the pit of kinetic energy substantially 
reduces the flux of clear air inside the cloud. 

Figure (c) shows the vertical derivative of the downward flux of clear air 
when Fr = 1.8 As the flow evolves, the downward flux reduces and its derivative, 
which represents the net variation of 1 — y at a given instant, rapidly tends to 
zero inside the cloud. This implies that the entrainment of clear air is conhned 
to a thin interfacial layer. 


4 Conclusion 

In this work we have carried out numerical simulations on the transport of en¬ 
ergy and scalars in a turbulent shearless mixing layer associated to temporal 
perturbation of the temperature lapse rate across the clear air - cloud inter¬ 
face. The perturbation locally introduces a stable stratification. This idealized 
configuration models some of the phenomena which are present in the kinetic 
dynamics of the cloud and clear air interaction, namely those linked to turbulent 
mixing and entrainment. 
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We have shown that this flow conflguration develops an horizontal layered 
structure characterized by a sublayer - with a kinetic energy lower than both 
the external regions - which acts as a barrier for the transport between the 
cloud and the external ambient. In our transient simulations, this flow structure 
appears when the buoyancy terms becomes of the same order of magnitude of 
the inertial ones, therefore the time needed for this transition becomes shorter 
when the stratification is more intense. Once buoyancy dominates and the 
new flow regime is reached, we observed two highly intermittent regions with 
opposite kinetic energy gradients. As a direct consequence, the entrainment is 
damped. 

Results obtained so far seemingly support the large eddy simulations of 
stratocumulus - topped planetary boundary layers. 
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